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Abstract 

in ' 

In this paper, a boundary integral method is used to solve an inverse linear heat 

f-H ■ conduction problem in two-dimensional bounded domain. An inverse problem of mea- 

Q_|. suring the heat flux from partial (on part of the boundary) dynamic boundary mea- 

^ fH surements is considered. An algorithm is given by using the fundamental solution, 

-i— > . 

a . 
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o\ : 

i> ■ 1 Introduction 

in ■ 

This paper is devoted to an inverse problem for a type of parabolic PDEs in a bounded 
two-dimensional domain. Here we consider initial boundary value problems for the 
heat equations by using a boundary integral approaches The inverse heat conduction 
problem arising in most thermal manufacturing processes has recently attracted much 
■ attention [21 [61 [5j [T71 [THl [20]. The typical case is the determination of the heat 

flux on an inaccessible boundary through measurements on an accessible boundary 
or inside the domain. Similar problems to ours have been studied by many authors, 
for example in one dimensional space we can refer to an approximate inverse method 
[9j [10] , a boundary element method [TBJ , a fundamental solution method j3[ [8] and 
some other method |19j . In this paper, we use the boundary integral method to solve 
this problem. This method uses the prescribed initial and boundary data, together 
with the fundamental solution of a given differential equation defined in some bounded 
domain, and we construct integral equations on the boundary of . In our case, the 
solution to the integral equation is a single layer potential. By the boundary integral 
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equation one can obtain the unknown kernel, and the solution to the given problem will 
be obtained by integrating the product of the fundamental solution and the unknown 
kernel over the boundary. The advantage of our approach is that the computation can 
be limited to the boundary, which reduces the problem from two dimensions to one 
dimension. As a result of the reduction, we may expect substantial savings in computer 
time and memory. The work outlined below is based on the use of single layer poten- 
tials. Ammari and Kang used boundary integral method to solve inverse conductivity 
problem and related problems [ilj. In [1], both the single layer potential and the double 
layer potential are used. In this article, our boundary integral method is based on the 
result of [14J, which gives a representation formula for the heat conduction problem 
with Neumann boundary condition. The equation is assumed to be homogeneous. The 
outline of our paper is as follows:... 

2 Problem formulation 

Let C R 2 be a bounded domain of M 2 . We denote dfl the boundary set assumed to 
be of class C . We denote by v the outward unit normal to fl on dfl. Let T > 0. We 
consider the following homogeneous heat equation: 

d t u - Au = 0, (x, t) G fi x [0, T] 

u(x, 0) = Uq(x), X £ ft , . 

u(x,t) = g(x,t), (x, t) G dfi x [0, T] ^' 

^(x,t) = <f>(x,t), (x,t) e on x [0,T] 

where u(x, t) is the temperature function, uq(x) the initial data, g(x, t) is a suitably 

du 

prescribed function, (j)(x, t) is the unknown heat flux and dtu = — is the rate of change 

at 

du 

of temperature at a point over time. Notice that d v u = — = Vu • v. 

ov 

Let 7(C) : [0, 1] — » M 2 be an analytic, 1-periodic function. We assume that the 
tangential derivative has the positive length |7'(C)| > for all < £ < 1. 

Throuout this paper we suppose that the closed smooth curve dfl in M? is param- 
eterized by the function j(C) as follows : 

an = {x = j(0, Ce[o,i]}. (2) 

The present paper proposes a boundary integral method for the numerical solution 
of the two-dimensional problem defined by ([1]). Our problem can be stated as follows: 

The inverse problem 

Let r C dtt denote a measurable smooth connected part of the boundary dfl (T = dfl 

du 

or not). The aim of this paper is to determine the heat flux (j)(x,t) = — — \r(x, t) from 

du 

measurements of : 

u(x,t) onTx[0,T]. 

For this purpose, we develop a boundary integral method .... as will be described in 
the next section. 
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3 Boundary integral method 

In this section we consider the basic boundary integral approaches for the solution of the 
initial boundary value problem ([1]). As in the time- independent case there are two main 
types of approaches, namely the direct method based on the representation coming from 
Green's formula and the indirect or layer methods (see for example [?]). We begin with 
the fundamental solution of the heat equation. In several spatial variables, the Green's 
function is a solution of the initial value problem (see for example [8] ) : 

d t G - AG = 0, and G(x, t = 0) = 5{x) (3) 

where 5 is the Dirac delta function. The solution to this problem in W 1 (n > 1) is the 
fundamental solution : 

\ Hit) . Ixl 2 , 

G(x,t) = T Z 7T7exp(- L - L ), 

v ; (4vrt)™/ 2 v At ' 

where H (t) is the Heaviside function and |x| = (x\ + x\ + • • • + x^) 1 ! 2 for x = 
(x 1 ,x 2 , ••• ,x n ) e W 1 . 

Let S and T> denote the classical single-layer and double-layer heat potentials : 

(Sq)(x,t) = / / G(x-y,t- s)q(y,s) da(y)ds, 
Jo Jan 

(Vip)(x,t) = / d v r y \G{x - y,t - s)<p(y,s) da(y)ds, 
Jo Jan 

and similarly we define V as the spatial adjoint of the double-layer T> and 7i the hyper 
singular heat operators as follows: 

(V'q)(x,t)= ! [ d u(x) G(x-y,t-s)q(y,s) da{y)ds, for (x, t) e dQ x [0, T], 
Jo Jan 

and 

{Htp){x,t) = ~ f ( d2G{ *~ V /~ S \ {y,s) da(y)ds, for (x,t) G 90x [0,T]; 



o J an du{x)d v ( y ) 
where q and (p sufficiently smooth functions. 

In terms of the operators S and T>, we introduce the following heat potential 

u = Sq-V(p inOx[0,T], (4) 

where u is the solution of ([1]). Now, by using the well known boundary behavior 
properties of single-layer and the double-layer heat potential the following result holds. 



Proposition 3.1 Let u be the solution of the problem Then, in both cases ho- 

mogeneous or inhomogeneous initial data uq, the Cauchy data ( , u|,9nx[o,T]> ^ u box[o,Tl) 
satisfy : 

u = Sq+(^I-V)v indnx[0,T] (5) 

8 u u= (^I + V')q + H(p mfflx[0,T], (6) 
where the functions q and <p are given as in Q). 
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Proof. For the case of inhomogeneous initial data uq = one can get the above results 
by using the well known boundary behavior properties of single-layer and the double- 
layer heat potential together with their normal derivatives. But, for inhomogeneous 
initial data uq ^ 0, the situation needs more explanations. It can be seen that a 
function 



satisfies 
and 



v(x,t) = / G(x - y,t)u (y) dy 
Jo Jan 



v t - Av = 



lim v (x, t) = uq(x). 



Setting w = u — v, we get 



d t w -Aw = 0, (x, t) G n x [0, T] 

w(x,o) =0, x g 

w(x, t) = g(x, t) - v(x, t), (x, t) G dil x [0, T], 

^( x ,t) = ( f ) (x,t)-^\(x,t), (x,t) e an x [0,T\. 



(7) 



Since w is the solution of ([T]) with homogeneous initial data, it can be solved to find 
the integral relations (JSD- dSl) • Thus we can solve the problem and get our results by 
the superpositions of v and the solution of homogenous initial data. □ 
Now, we introduce the following notations and the anisotropic Sobolev space to be 
used in the squeal. A rather comprehensive of the basic presentation of these Sobolev 
spaces the reader can see [HI ?]. For given r, p > 0, we have the space 

H r ' p (Q x [0,T]) = L 2 ([0,T];F(!]))nff p ([0,T];L 2 (O)). 

The space H r,p (dQ x [0, T]) is defined analogously by replacing f2 by d£l. Moreover 
the following subspace of H r ' p (d{l x [0, T]) is well defined 

ff' p (dn x [0,2*1) = {v = w\ dnx[0 , T] : w G H r ' p (n x M), w(-,t) = 0,t < 0}. 

The norm of the space H r ' p (Cl x [0,T]) is denoted by || • ||r,p ; nx[o,T] an d the norm of 
h r > p (dn x [0,2*]) is denoted by || • \\ r , P ;dnxlo,T\- 
Next, introducing the useful spaces 

£ = {v G H 2 ^(nx[0,T}) : (A + d t )v G L 2 (n x [0, T}), £U| anx[0jT] = 0,v(-,T)\ n = 0}, 

X = h 1/2 ' 1/4 (dn x [0,T]) and the associate dual space X' = fT l/2 " 1/4 (<90 x [0,T]). 
Then the following definition appears. 

Definition 3.1 For given f G we say that u is a weak solution of (T7]j if u G 
2f 1,:L / 2 (r2 x [0,T]) and satisfies the following duality product: 

(u,Aip + d t ip) = -(f,ip)dnx[o,T], 

for ijj G £. 

In terms of the last notations, the following mapping properties of the single- and 
double-layer heat operators holds. 
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Lemma 3.1 Let the operators S and T> be defined as in section 2. Then the fallowings 
hold: 

(1) The single-layer heat operator S : ff'^idtt x [0,T]) -> tf+ 1 M 2 ir+ 1 ) (<9fi x [0,T]) 
is an isomorphism for all r > —1/2. 

(^J T/ie operator \L + D : X' ^ X' is an isomorphism. 

Proof. The claim (1) can be inspired directly from Theorem 4.3 in [?]. The claim (2) 
follows by a little modification from [TJ. □ 
Now we're ready to prove the following result. 

Theorem 3.1 Assume that (fi G X', Then the function u G H 1 ' 1 / 2 ^ x [0, T}) is a weak 
solution of Ul) if and only if u has the representation ^ such that ip G X' solves the 
equations f5p-(0|). 

Proof. Let (fi G X' , and let p G X' be the unique solution of ©-(JS]). Then, by 
combining relation ([5]) and the direct representation we may get the following 
boundary integral equation of the second kind : 

(ll + V)p = S(f>. (8) 

Since the set D(<9$7 x [0, T]) is dense in the space X' , we can choose a sequence ip n G 
D(<9$7 x [0,T]) such that p n —> tp in X', and so the following sequence (fi n = + 
T>)cp n is well defined. By the mapping properties of S and of the operator ^1 +T> found 
in Lemma |3.1| the function <j) n is also a smooth function of dVt x [0, T]; moreover we 
have (fi n — > 4> in X' . 

Now, let u n be the corresponding classical potential 

Un = Sip n - Vcfi n , 

which by the construction satisfies 

dvU n \dny.[$,T] = (fin- 

Obviously, u n is a weak solution of (pQ) with the Neumann data <fi n . On the other hand, 
As done for the Dirichlet-type initial boundary value problem in [?], we can conclude for 
our problem that the mapping (fi — > u, is continuous and we have u G H 1 ' 1 / 2 ^ x [0, T]) 
such that 

IMIi,i/2,nx[0,T] < c||0||i/ 2 ,i/4,anx[O,T], where c is a positive constant. (9) 

Hence, by the continuity Q we have the convergence u n — > u in ff 1,1 / 2 (J7 x [0,T]). 
Next, we define tfi = Sip - Vcfi in Q, x [0, T] , then for all v G D(0 x [0, T]) it follows that 

(u,v) = lim (-u n ,f) = lim ((Sp n ,v) - (D(fi n ,v)) = (ifi,v), 
which implies that u = ip = Sip — Vcfi. □ 
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4 Numerical scheme for the inverse problem 

In this section we propose a numerical method to solve our inverse problem. The 
numerical method is based on the boundary integral equation in Proposition ??. As the 
measured data for inverse problem, the numerical data obtained by solving the direct 
problem can be used. The inverse problem is then to solve the following problem: 

d t u-Au = 0, (x, t)£flx [0, T] 

u(x,0)=0, x£fl 

u(x,t) = g(x,t), (x, t) e dQ x [0, T] { ' 

^(x,t) = <f>(x,t), (x, t) E dfl x [0, T]. 

We use the following equation given by (J3|) 

u \anx[o,T\ =Sq-Vtp inOx[0,T], (11) 

together with the relations ©-([6]) to solve the above problem. 

To give our numerical numerical method, we introduce the following result. 

Lemma 4.1 Let <p = [u]qq x [q : t], u solution of MO)) . In term of the hyper- singular heat 
operator the heat flux is solution of : 

Hip = <P, (12) 
where [w]anx[o,T] means the jump of the function u via the boundary d£l x [0, T\. 

The boundary integral equation of the first kind given in Lemma 14.11 is deduced from 
a normal derivative applied to the double layer representation which itself given by 
managing relation (JH) into relations given in Proposition 13. 11 

To proceed with our numerical scheme, we may follows two cases, measure on the 
boundary d£l and measure on a smooth connected subset of the boundary. 



4.1 Measure on the boundary dQ 

The aim of this section is to reconstruct <j)(x, t) from measurements of u(x, t) on the 
boundary dfl x (0, T). For this purpose, we develop the subdivision of [0, 1] : 

Co = 0, d = Co + ih, i = 1, 2, • • • , N, where N is an integer and h = l/N. (13) 

Analogously, we assume the subdivision of [0, T] : 

t = 0, tj = t + jti, j = 1, 2, ■ • ■ ,N' where N' is an integer and h' = T/N'. (14) 

Then, by using Lemma 14. II the following main result follows. 

Theorem 4.1 Let g G X be a given function and u € if 1 ' 1 ' 2 ^ x [0, T]) be the solution 
of fTOj) . Assume that we have the subdivisions [T3\) -p4\). Suppose that the heat flux of 
the problem ilO\) is continuous up to the inner side of the inaccessible boundary (90, 
then the unknown data <f> can be recovered by the following discrete scheme 

ma*) = - g g (t ._ t0 2 mCt)l • [ - ^ (oo+ (is) 
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>(V(« • (7(&) - 7(o))) (7(c ^:; (a)) ] ex P( - '^)-_^>i 2 ), 

for i = 1, 2, • • • , N; j = 1, 2, • • • , N' and 7 is given by Section 2. 

Proof. Let u be the solution of the inverse heat problem (|l(jp and inserting the 
expression of the hyper-singular operator 7i into relation (|12p . we get 

</>(x,t) = ~ f I &2G } X ( ~J 7t ( ~, S) 9(y,s) da(y)ds, for (x,t) G dSl x [0,T]. 
Jo Jdn dv{x)dv(y) 



By change of variable x = 7(C) 5 C £ [0) 1]) we write 
'0 Jo du(x)du(y) 



0(7(0,*) = -/ / a , ^ , s (7(C)-7(C0, *- s )g(y> g)IY(C0l rfC^s, for (C,t)G[0 5 l]x[0,T]. 



(16) 

Then according to (|13p -(|14p. we can discretize relation (|16p as follows : 



0(7(0), ^) = -fefe'EE ^ ^? ) (-y(Ci) -7(a),^ -*i)g(7(Cfc),*oiV(Cfe)i- (17) 
fe=i ;=i ^ ^ 

On the other hand the normal derivative of the fundamental solution G of ([3]) in two 
dimensional space is : 

d u{y) G{x-y,t-s) = 4(t _ g)2 ex P(- (t _ g J - 

Thus, we can derive a gain by v{x) to get throw the formula V(A ■ B) = A x (V x 
B) + B x (V x A) + (A- V)B + (B • V)A that : 

<9 2 G 1 
du(x)du(y) iX -y^- S ^ = ~Mt^W V ^ • & " X) X (V X U{V)) 

HHV) ■ V)(y - x) + ((y - x) • V)«/(lO + 2(«/(l0 • (y - x)) exp(-^^). 

' I — S J (t — 8) 

Therefore to achieve the proof, we insert the last formula into relation (|17|) by taking 
x = 7(0) and y = j{Ck)- □ 



Now, one can give an approximation to the solution u(x, t) of the problem (jlOp 
by inserting the discredited formula of the heat flux given by Theorem 14.11 into the 
representation as : 



N N' 



i( Xi ,tj) :=^ , ^^G( 7 (O)-7(a),*i-^)0(7(Cfe),^)l7 / (a)| (18) 



k=l 1=1 

N N' 



Fig 

-^TE^(^ i )-7(a) ! ^-^(7(a),^)i7'(a)i 
k=i 1=1 
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4.2 Measure on a smooth subset of the boundary 

Let Tec d£l denote a measurable smooth connected part of the boundary d£l and T c 
denotes dCl \ T. Introduce the trace space 

X =|u G X,v = on T c x (0,T)|. 

Here and in the sequel we identify g defined only on T with its extension by to all 
dfl. 

The aim of this section is then to identify the heat flux <f> from the local measure on 
the Cauchy data g G X . To do this, we may assume that there exists C* £ (0, 1) such 
that 

r := {7(0 : C € [0, C*], where C* « 1}- 
As done in last section, we introduce the subdivision of [0, C*] : 

Co = 0, Qi = Co + ifj i = 1) 2, • • • , M where M is an integer and r = C*/M. (19) 
Similarly for [0, T], 

t = 0, tj = t + jr', j = 1, 2, • • • , M' where M' is an integer and r' = T/M'. (20) 
As done in Theorem 14.11 for the local measurement we have the main result. 

Corollary 4.1 Let g G X be a given function and u G H l,1 ^ 2 (Q, x [0, T\) be the solution 
°f Iil0\) - Assume that we have the subdivision Uty -Hi ty) . Suppose that the heat flux 
satisfy the hypothesis in Theorem \4-l\ Then the data 4> c o- n be measured on T as : 

mCi),tj) = -rr'J^T, du( )dv( M d) -7(Cfc),ti -tiMl((k),ti)W((k)\ (21) 

h=l 1=1 ^J> 
T(* v^y^ g(7(Cfc)»^) l'(Ci) r , (t x. 
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2(7'(C,) • (7(C) - 7(C,))) (7( ^:; (Cfc)) ] ex P ( J^7^i 2 ), 

/or i = 1,2,--- , M; j = 1, 2, ■ • • , M' and 7 is given by Section 2. 

Proof. Consider that u is the solution of the inverse heat problem (|10p . then by (|12p . 
we can write 

<f>{x,t) = -fl &2G } X -y> t - s) g(y,s)do-(y)ds, for (*, i) G T x [0, T]. 
Jo Jrur c dv(x)dv{y) 

The fact that <? G Af we can reduce 

<P(x,t) = - f f d2 °} X ~ V /~ s) g{y,s) da(y)ds, for (x,t) G T x [0,T]. 
Jo Jr du(x)du{y) 



By change of variables and by considering the subdivisions (|T9]) - (pO|) one can deduce 
from the proof of Theorem 14.11 that 

M M> 

matj) = -rr^^ (7(0) - 7(a), ^ - tiMimtiMm. 

fe=l 1=1 ^ ' 

Then, the proof achieves by inserting in last relation the possible normal derivations 
of the function G. □ 
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4.3 Numerical examples 



Numerical evaluations of the heat flux 4>(x, t) are now obtained by solving the integral 
equations (jlip and (|12p with the specific parameterization and discretisation of the 
boundary. 

In this section, we assume that the curve is parameterized by 7(C) = (cos(2-7r£), sin(27r£)) 
and we suppose that T = 10, C* = 10~ 2 . 

4.3.1 Measure on dtt x [0,T] 

We consider the following direct problem : 

d t u-Au = 0, (x, t) E Q x (0, T) 

u(x,0)=0, xen , s 

u(x,t) = 2|ic|cos(3i) J (x, t) € <9S7 x (0, T) 1 ' 

^(x,t) = ( f ) (x,t), ( x ,t)ednx(o,T). 

We solve this problem by using Theorem 14,11 given by boundary integral method which 
is presented in sections 3 and 4. The parameter is chosen as N = 50, N' = 100. A 
numerical result is shown in Figures 

4.3.2 Measure on T x [0, T] 

For the case of smooth subset T we consider the following problem : 

d t u-Au = 0, (x, t) E Q x (0, T) 
u(x,0)=0, ieO 
u(x,t) = 2\x\cos(3t), (x, t) e T X (0, T) 1 ' 

^(x,t) = 4>(x,t), (x,t) erx (o,T). 

This problem can be solved by using Corollary 14.11 The parameter is chosen as 
M = 50, M' = 100. A numerical result is shown in Figures 

4.4 Conclusion 

A boundary integral method for the two-dimensional inverse heat conduction problem 
is considerably discussed. We presented a numerical scheme for the inverse problem. 
The heat flux was measured from the whole boundary and from a smooth subset of this 
boundary. In this paper we restricted our selves to homogeneous conduction problem, 
but the case with an external source may be considered in a forthcoming work. 
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